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Abstract 

In order to protect the equipments and make rules for 
lightning protections, it's essential to record accurate 
lightning overvoltage waves. The use of serial reactors which 
are intended to limit large current brings new problems for 
measuring the lightning overvoltage. This paper analyzes 
the change of the lightning overvoltage through the serial 
reactor and the bus, and concludes that the overvoltage at 
the entrance of the transformer can't reflect the invasion 
voltage of the substation well. In review of the position of 
the arrester in the substation, a method for calculating the 
lightning overvoltage is presented by using the leakage 
current of the arrester, and reverse iterative method is put 
forward in the calculation of the PINCETI model. The 
leakage current of HOkV arrester is simulated with ATP- 
EMTP and the amplitude and steepness of the lightning 
overvoltage are calculated, which verified the effectiveness 
of the proposed method. 
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Introduction 

With the development of power system, the height of 
the overhead transmission lines continues to increase 
and the possibility that power outages are caused by 
lightning is also increasing. Accurate lightning 
overvoltage waveforms can help us to analyze the 
reason of accidents, and improve insulation 
coordination [1-2]. 

There are a few ways to get invaded overvoltage 
waveforms, but they have some defects. PTs and CVTs 
don't have excellent transient characteristics [3], 
capacitance dividers may bring potential dangers [4-5], 
extra voltage sensors may cause the ground line break 
[6], and the method that the voltage is calculated 
according to leakage current doesn't reflect the 
invaded voltage [7]. 

The use of serial reactors can limit short circuit 
currents, while they change the propagation 



characteristics. This paper analyzes the influence of 
bus and serial reactors and simulates it using ATP- 
EMTP. 

MOA (Metal Oxide Arrester) are widely used in 
power system to limit overvoltage and protect 
electrical equipments for its excellent performance. 
Because the arrester is the first equipment that the 
overvoltage goes through, a method to calculate 
lightning overvoltage using the leakage current is 
presented and is verified by ATP-EMTP simulation 
results. 

The Propagation of Overvoltage In Station 

The capacitive and inductive equipments, the 
distributed capacitance of equipments and the mutual 
capacitance between equipments can impact on the 
spread of the lightning overvoltage [8-9]. 

1) The Effect of Capacitance 
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(A) THE SCHEMATIC DIAGRAM OF LIGHTNING 
OVERVOLTAGE THROUGH CAPACITANCE 
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(b) The equivalent circuit 

FIG.l PROMULGATION CHARACTERISTICS OF LIGHTNING 
OVERVOLTAGE WAVE AT THE SHUNT CAPACITOR 

As shown in Fig.l, when the lightning overvoltage 
passes the capacitor C, according to Peterson Law, 
the refracted voltage can be calculated from 
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formula (1): 



U 2q ~ hq^l ~ 
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Where: T = Z,Z 2 C I {Z x +Z 2 ) . 
2) The Effect of Inductance 
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(A) THE SCHEMATIC DIAGRAM OF LIGHTNING 
OVERVOLTAGE THROUGH INDUCTANCE 
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(B) THE EQUIVALENT CIRCUIT 

FIG.2 PROMULGATION CHARACTERISTICS OF LIGHTNING 
OVERVOLTAGE WAVE THROUGH SERIAL INDUCTOR 

When the lightning overvoltage passes the serial 
inductor L, the refracted voltage can be calculated 
according to formula (2) 
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Where: T = L/(Z l +Z 2 ). 
3) ATP-EMTP Simulation 

A model of HOkV substation is built as Fig. 3 [10]. It 
has three lines and Line3 has serial reactor L. 
Where: the reactor L=3H and the distributed 
capacitance of the bus C=0.006 u F. Line 3 is struck 
at 0.001s and the voltage waveforms monitored are 
presented in Fig.4. (ua is the voltage at arrester and 
ut is at the bolster of the transformer). 
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FIG.4 THE LIGHTNING OVERVOLTAGE WAVE AT THE 
POSITION OF THE ARRESTER AND THE TRANSFORMER 

In Fig.4, the voltage waveforms are very different 
in amplitude and steepness monitored at the two 
places and the voltage at the bolster can't reflect the 
invaded voltage of the substation. So, a method to 
calculate the characteristic parameters using 
leakage current of MOA is posted. 

Calculation of the Characteristic Parame- 
ters of Lightning Overvoltage 

The amplitude and steepness are the most important 
parameters of lightning overvoltage, which can direct 
to analyze the cause of malfunction and make 
protection rules. The leakage current can be measured 
by current sensors and a method named reverse 
iterative method is posted to calculate overvoltage. 

1) The Measurement of Leakage Current 

As shown in Fig. 5, the system includes three parts 
[11-15]: current sensor, data acquisition and 
transmission module and data analysis module. 
When the voltage across arrester is too high, the 
leakage current gets high and the relationship 
between voltage e measured by secondary coil 
induced in coil and the leakage current can be 
calculated as formula (3)-(4) 



_ dy/ _ jy^dB _ NSju Q ^ di 
dt dt Inr dt 



2nr 



jedt 



(3) 
(4) 



Where: 



FIG.3 THE SCHEME OF A 500KV SUBSTATION 



e — voltage induced in the secondary coil; 

i — leakage current; 

N — turns of secondary coil; 

B — flux density around the ground line; 

S — area of secondary coil; 

r — radius of secondary coil; 
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V — flux in secondary coil; front current surge (rl/T2 u s); Ur8/20 represents 

the residual voltage at lOkA fast front current surge 

Mo ... permeability of air. with a r8/20 U s shape; Un is the rated voltage of the 

arrester. 
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FIG.5 ON-LINE MONITORING SYSTEM OF LEAKAGE CURRENT 
OF ARRESTER 

2) MOA Models 

Two models of MOA are shown in Fig.6, IEEE 
model and PINCETI model. The parameters of 
IEEE model are more difficult to determine, and 
the parameters of PINCETI model only depend on 
test parameters [16-18]. In this paper, the 
overvoltage across the arrester is calculated on the 
base of PINCETI model. 



(a) IEEE Model 



R u A 



(b) PINCET Model 
FIG.6 SIMPLIFIED MODEL OF ZINC OXIDE ARRESTERS 

In PINCETI model, R equals 1MD; The u~i 
characteristics of nonlinear resistances AO and Al 
are shown in Fig. 7; To calculate LO and LI, 
equations (5)-(6) can be used: 



■1/T2 ^r8/20 ) 

^=(y n Mu ryT2 -u rs/m )/(ij rwx )jj n 



(5) 



(6) 



FIG.7 THE CHARACTERISTICS OF AO AND Al (THE VOLTAGE IS 
IN P.U. REFERRED TO THE VR8/20) 

3) The Mathematical Model of the Arrester 

In PINCETI model, equations (7)-(ll) can reflect 
the relationship between the voltage across the 
arrester and the leakage current. 
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(11) 



Where: fO > fl represent the u~i function of 
nonlinear resistances AO and Al; u is the wanted 
variable. 

4) The Calculation of Voltage 

A method named reverse iterative method is put 
forward in order to calculate the voltage. The basic 
idea of reverse iteration method is: in Fig.8, if x and 
y has the same variation, y has only one value in 
correspond with x. According to the errors between 
calculated and measured data, a penalty factor is 
introduced and y can be recalculated until the error 
is within the allowable range. 



Url/T2 represents the residual voltage at lOkA fast 
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FIG.8 THE CIRCUIT CONTAINING NONLINEAR RESISTANCE 
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In Fig. 9, f represents the function between the input 
and output data. If f is monotonous, x has one 
value to satisfy that y equals yref [19]. 
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FIG.9 ANALYSIS OF BACKWARDS ITERATIVE METHOD 

In Fig.6 (b), if il increases, i also increases. The 
parameters in PINCETI model can be calculated 
using reverse iterative method. 

At any time k, the specific solution steps are as 
follows. 

CD Divide u~i curves of AO and Al into n segments, 
discrete formulas (9)~(11) and assume that each 
segment is linear. 

© Assign il an initial value il(k) and determine 
which segment il(k) belongs to. If il(k) belongs to 
ml, then the voltage across the resistance Al is 



"i (k) = fi (h (k)) = u lml0 + (/, (k) - 

! 'lmlo) * ("lrnll _ "lmio) I (hmll ~ ! lmlO ) 



(12) 



Where: (ilmlO, ulmlO) , (ilmll, ulmll ) 
represent the start and the end of segment ml. 

® Calculate voltages across the nonlinear resistor 
AO and the inductor LI : 



«Ll(*) = A' 



At 



u (k) = Mj(fe) + u Ll (k) 



(13) 
(14) 



i If uO(k) belongs to mO, the current through AO is 

i (k) = i OmO0 + (u Q (k) - u 0m00 ) • 

(15) 



Omul 'OmOO 



)/(«0 



Where: (iOmOO, uOmOO) , (iOmOl, uOmOl ) 
represent the start and the end of segment mO. 

© Calculate the current through LO and the voltage 
across LO : 



io( k )-'o( k - 1 ) 



U Lo( k ^ ~ A) " 



At 



(16) 
(17) 



u(k) = u L0 (k) + u Q (k) 

© Calculate leakage current 

i(k) = i (k) + u(k) / R 



(18) 



(19) 



® If the error between i(k) and iref (k) is not within 
the allowable range, according to equation (20), 
correct il(k), and repeat the above calculation 
processes until the error is small. Where a is the 
penalty factor. 

i, (k) = i, (k) + a» (i(k) - i ref (k)) 

The value u in process © is what we want. 

EMTP Simulation 

1) Parameters of Arrester 

A HOkV transmission system is simulated using 
ATP-EMTP. The type of MOA is Y10W5-100/248, 
and its electrical characteristics are shown in Tab.l 
[20]. (voltages are in kV) 

Tab.1 Measured electrical parameters of MOA: Y10W5- 
100/248 
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In Tab.l, Un is the rated voltage of the system; Ur 
is the rated voltage of MOA; Uc is the continuous 
operating voltage of MOA; Url/T2,Iamp represents 
the residual voltage for surge current impulse with 
lamp amplitude and rl/T2ps shape. 

The characteristics of AO and Al in PINCETI model 
are in Tab. 2. 

Tab.2 V-A characteristics of AO Al 



I/kA 
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i Calculate lightning voltage u 



2) Results 

After the line was struck by lightning, the arrester 
discharged to release lightning energy. And the 
leakage current measured is shown in Fig.10. 
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FIG.10 THE LEAKAGE CURRENT OF THE ARRESTER 

In Fig. 11, the voltages calculated and measured 
across the arrester in shown (blue dotted line 
represents the calculated data and the red line 
represents the measured data). 
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FIG.ll THE COMPARISON BETWEEN THE CALCULATED DATA 
AND THE TRUE DATA 

The two parameters amplitude and steepness of 
the lightning overvoltage are in Tab. 3. The 
steepness is represented by the time of rise. 

Tab.3 The comparison of the characteristics between the 
calculated data and the true 

Calculated Measured Error 

data data 
Amplitud 2.2676*10 5 V 2.2733*10 5 V 0.25% 
e 

Time of 3.505xl0-"s 3.503*10- 4 s 0.057 
rise % 

In Tab.3, the parameters of the lightning 

overvoltage can be calculated correctly with the 

method presented in this paper. 

Conclusion 

This paper simulates and analyzes the impact of serial 
reactor and bus, and concludes that the voltage at 
bolter of transformers can't reflect the lightning 
overvoltage invaded; In order to calculate the 
lightning overvoltage, a method named reverse 
iteration method is posted; The voltage calculated 
with this method is very similar to the real voltage, 
and this method is verified effective. 

The method presented in this paper can calculate the 



amplitude and steepness of the lightning overvoltage 
correctly, which is helpful to analyze malfunction 
reason and make protection measures. 
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